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ABSTRACT 

We present a model to explain the wide range of abundances for heavy r-process elements (mass 
number A > 130) at low [Fe/H]. This model requires rapid star formation and/or an initial population 
of supermassive stars in the earliest condensed clots of matter to provide a prompt or initial Fe inventory. 
Subsequent Fe and r-process enrichment was provided by two types of supernovae: one producing heavy 
r-elements with no Fe on a rather short timescale and the other producing light r-elements [A < 130) 
with Fe on a much longer timescale. 

Subject headings: Galaxy: evolution — Stars: abundances — Stars: Population II 



1. INTRODUCTION 

We present a phenomenological model for the abun- 

' dances of Fe and heavy (mass number A > 130) and 

I light {A < 130) r-process elements (r-elements) in very 

■ metal-poor stars. These stars formed early in Galactic 
' history when only a small number of massive stars had 
. evolved to become Type II supernovae and added heavy 
' elements to the interstellar medium (ISM). Observations 
; (e.g., McWiUiam et al. 1995; McWiUiam 1998) show that 
. there is a lack of correlation between the "metallicity" 
' measured by [Fe/H] and abundances of r-elements above 
' A ^ 135 (e.g., Ba and Eu) for stars with -3.1 < [Fe/H] 

• —2.5. This disagrees with the view that Fe and r-elements 

' are always coproduced by supernovae. Our approach to ac- 

, count for this result is based on a two-component r-process 
model (Wasserburg, Busso, & GaUino 1996; Qian, Vogel, 

[ & Wasserburg 1998; Qian & Wasserburg 1999, hereafter 

. QW99) which attributes heavy r-nuclei to high frequency 

' supernovae (H events) and light r-nuclei to less frequent 

[ ones (L events). We propose that prompt production of 

■ Fe with minor coproduction of heavy r-elements in the 

■ earliest era of the Galaxy first enriched the ISM up to 
, [Fe/H] ~ —3. Subsequent production oi A > 130 nuclei 

with negligible coproduction of Fe by H events then re- 
sulted in a wide range of abundances for heavy r-elements 
at —3.1 <^ [Fe/H] ^ —2.5. We associate further Fe enrich- 
ment of the ISM for [Fe/H] > -3 with L events. Addition 
of H and L events over a sufficiently long timescale then led 
to a correlation between [Fe/H] and abundances of heavy 
r-elements at [Fe/H] ^ —2.5. 

Mathews, Bazan, & Cowan (1992) discussed using 
Galactic chemical evolution to constrain the site of the r- 
process. A number of recent studies (Ishimaru & Wanajo 
1999; Tsujimoto, Shigeyama, & Yoshii 1999; McWiUiam 
& Searle 1999; see also Raiteri et al. 1999) focused on the 
relationships between [Fe/H] and abundances of heavy r- 
elements in the early Galaxy. A common consensus is that 
chemical enrichment of the ISM at very early times was 
grossly inhomogeneous and diverse yields of individual su- 
pernova events had dramatic effects on abundances in very 
metal-poor stars. Our present work differs from previous 
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studies in that the diversity in r-process production by 
supernovae is introduced through the two-component r- 
process model based on solar system data independent of 
the stellar observations. Further, we propose a prompt 
mechanism for Fe production considering special condi- 
tions of the early Galaxy. In §2 we describe the framework 
of two r-process components. In §3 we use this together 
with a postulated prompt Fe source to explain the obser- 
vational results on abundances (especially for Ba and Eu) 
in the early Galaxy. We discuss prompt Fe production and 
give conclusions in §4. 

2. TWO i?-PROCESS COMPONENTS 

Observations by Sneden et al. (1996, 1998) demon- 
strated that abundances of r-elements in the Pt peak 
{A - 195) and down to Ba {A ~ 135) in CS 22892- 
052 ([Fe/H] = -3.1), HD 115444 ([Fe/H] = -2.77), and 
HD 126238 ([Fe/H] = -1.67) are in remarkable accord 
with the solar system r-process abundance pattern (the 
solar r-pattern). Assuming that a single r- pattern ex- 
tends from Ba to the actinides above the Pt peak, Cowan 
et al. (1997, 1999) discussed using the Th/Eu ratio 
as a Galactic chronometer. However, the discovery of 
^^^Hf (lifetime fi82 — 1.30 x 10^ yr) in meteorites with 
(•I82jjf/i80jjf)ggp = 2.4 X lO-'^ (Harper & Jacobsen 1996; 
Lee & Halliday 1995, 1997) at the time of solar sys- 
tem formation (SSF) provided a new twist to our un- 
derstanding of the r-process. Although both ^^^Hf and 
^^^I (fi29 = 2.27 X 10^ yr) are produced essentially only 
by the r-process, the ^^^Hf data and abundance ratio 
(I29l/i27i)ggp = 10-4 (Reynolds I960; see also Brazzle 
et al. 1999) cannot be explained by a single type of r- 
process events. Based on this, Wasserburg et al. (1996) 
concluded that there had to be at least two distinct types 
of r-process events: one (H) occurring on a timescale 
Afj ^ iO^ yr, commensurate with that for replenishment 
of a typical molecular cloud with fresh supernova debris, 
and the other (L) occurring on a much longer timescale 
Al ~ 10^ yr. They further pointed out that relative to 
the solar r-pattern, there should be frequent abundance 
excesses of r-elements in the Pt peak over those in the 
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A ^ 130 peak below Ba at low metallicities. Qian et al. 
(1998) showed that in a two-component model to account 
for the solar r-pattern, the total mass yield of an L event 
must be ~ 10 times that of an H event. They also showed 
that it is not readily possible to produce the A ~ 130 
peak without substantially populating the region beyond 
this peak. They further speculated that H events might 
be associated with production of black holes and L events 
with production of neutron stars. 

A preliminary report on the observed Ag {A 107) 
abundance in CS 22892-052 by Cowan & Sneden (1999) 
appeared to support the meteoritic prediction that more 
than one type of r-process events may be required. Based 
on this, Qian & Wasserburg (QW99) sought to establish a 
quantitative basis for predicting the yields of r-elements in 
H and L events. Their model uses the following assump- 
tions: (1) each H or L event has fixed r-process yields; (2) 
the products of individual H and L events are mixed with 
a "standard" mass of hydrogen in the ISM; and (3) the 
solar system inventory of all r-nuclei (both stable and ra- 
dioactive) are the result of production by H and L events 
over a time Tup ~ 10^" yr preceding SSF. The timescales 
for occurrence of H and L events in a volume correspond- 
ing to the standard mixing mass of hydrogen are Ah and 
Al. Under the above assumptions, at the time of SSF, the 
number of nuclei for a stable nuclide S (e.g., ^^^I or ^^^W) 
mixed with a standard mass of hydrogen is 



SSF 



H^UP 



YkT, (1) 



while that for a radioactive nuclide TZ with Typ 
(e.g., 1291 or i82Hf) is 



SSF) 



^fexpHMnMO ^ 

1 - exp(-AH/fK) exp(AL/fTC) - 1 



yL 



(2) 



In equations (|l|) and (||) or denotes the yield in an 
H or L event. The yield for a radioactive species is taken 
to be about the same as that for the corresponding sta- 
ble nuclide. In equation (|^) the interval between SSF and 
the last L event is assumed to be Al, while that between 
SSF and the last H event is specified by i? (0 < < 1). 
From the meteoritic data on ^^^l and ^^^jjf^ .^y^g found 
that Ah < 1.85 x 10^ yr and Al > 1.06 x 10^ yr for 
19 = 1 (scenario A). Further, the yields for ^-^"^I and ^^^W 
in an L event relative to those in an H event, Y^^/Y^^ 
and Y1^2/Yi82j rather strongly constrained for all cases 
bracketed by scenarios A and B {i) = 0). Given t?. Ah, and 
Al, the yields for ^^^j ^nd ^^'^W in an H or L event can 
be obtained from the relevant solar system data through 
equations (|l|) and (||). 

To generalize the results to other r-nuclei, Qian & 
Wasserburg (QW99) assumed that the yield template for 
A > 130 nuclei associated with ^^'^W is the same as the cor- 
responding solar r-pattern for both H and L events based 
on the observations by Sneden et al. (1996, 1998). The 
yields for A < 130 nuclei associated with i27j jj^ q^^i H 
or L event were also chosen to follow the corresponding 
solar r-pattern. Consequently, the yield for a stable nu- 
cleus S with mass number A in an L event relative to 
that in an H event is Y^/Y^ = Y^^j/Y^^ for A < 130 or 
yL/yH = Y]ig2/^i82 for A > 130. The fraction of S nuclei 



contributed by H or L events to the corresponding solar r- 
process abundance are -F;" (5) = l/[l-t-(l^/Y|^)(AH/AL)] 
or F^{S) = 1 — F^{S). From equation (|l|) the abundance 
of stable S nuclei resulting from a single H or L event 
contaminating a standard mass of hydrogen in the ISM 
is logeH(5) = log €0,^(5) + logF^iS) - log(Tup/AH) or 
logeL(5) = logeo,,(5) + logi^^^^) _ iog(Tup/AL). Here 
the spectroscopic notation log e(5) = log(iS/H)-|-12 is used 
{S/Ti. being the number abundance ratio of S to hydro- 
gen). Thus, given i?. Ah, and A^, there is a quantitative 
prediction for the abundances resulting from a single r- 
process event. For example, a single H event gives rise to 
logeH(Eu) « -3.0 to -2.2 and logeH(Ba) ^ -2.1 to -1.3 
over a wide range of model parameters. The results from 
a mixture of multiple events can be calculated simply by 
adding the number of nuclei produced in each event and 
then converting it to the corresponding log e value for the 
mixture (see QW99). 

3. THE IRON CONUNDRUM 

Observational data (Gratton & Sneden 1994; 
McWiUiam et al. 1995; McWiUiam 1998; Sneden et 
al. 1996, 1998) on Ba/Eu, loge(Eu), and [Fe/H] for 
low-metallicity stars are shown in Figure 1. As can 
be seen, there is a wide dispersion in log e(Eu) at 
—3.1 <^ [Fe/H] ^ —2.5 while Ba/Eu is essentially con- 
stant. The r-process accounts for over 90% of the solar Eu 
inventory but only about 20% of the solar Ba inventory 
(Kappeler, Beer, & Wisshak 1989; Arlandini et al. 1999). 
Consequently, the clustering of Ba/Eu around the solar 
r-process value exhibited in Figure 1 confirms the earlier 
proposal by Truran (1981) that heavy elements such as 
Ba in very metal-poor stars originated from the r-process. 
There are few Eu data at —4 <^ [Fe/H] ^ —3. However, 
sufficient Ba data at these metallicities (McWiUiam et al. 
1995; McWiUiam 1998) are available and shown in Figure 
2 (region A) . Both Ba and Eu data show a wide dispersion 
at -3.1 ^ [Fe/H] -2.5 (region B) and suggest that Fe 
and heavy r-nuclei arc not coproduced by common super- 
novae (H events). The loge(Ba) and loge(Eu) values from 
a single H event in our model (QW99) arc indicated by the 
zone marked "1 H" in the corresponding figure. We expect 
that increases of loge(Ba) and loge(Eu) above this zone 
are dominantly the results of adding more H events to the 
standard mixing mass of hydrogen. For example, if we 
take loge(Eu) « -2.5 for HD 122563 ([Fe/H] w -2.7) as 
representative of a single H event, then the Eu abundance 
in CS 22892-052 ([Fe/H] « -3.1) would correspond to 
- 30 H events. As the Fe abundance in CS 22892-052 is 
smaller than that in HD 122563, it is evident that Fe can- 
not be significantly produced by H events. Thus the wide 
dispersion in Ba and Eu abundances at [Fe/H] ^ —3 poses 
a conundrum of Fe production in the early Galaxy and 
suggests that [Fe/H] is neither related to heavy r-element 
production nor a reliable chronometer (see QW99). 

We find that the Fe conundrum can be resolved by pos- 
tulating an initial or promptly-generated Fe inventory that 
existed before the occurrence of H and L events. By "ini- 
tial" we mean very early stages during which an inven- 
tory of Fe was provided to the ISM from which regu- 
lar stars later formed with no other temporal connection. 
Prompt Fe production is considered to be associated with 
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cocvolution of all stars from an initial gas clot with no 
Fc. In both cases the mechanism for Fe production must 
have ceased at [Fe/H] ~ —3. Then non-Fe-producing H 
events and Fe-producing L events (see below) began to 
occur. The frequent occurrence of H events would result 
in a range of abundances for heavy r-elements such as 
Ba and Eu at [Fe/H] ~ —3 while the less frequent oc- 
currence of L events would lead to increases in Fe abun- 
dance above [Fe/H] ~ —3. A correlation between [Fe/H] 
and abundances of heavy r-elements would then be es- 
tablished when sufficient Fe was produced by L events to 
overwhelm the inventory produced by the initial/prompt 
Fe source which was not related to "typical" supernovae. 
The existence of some stars with —4 <^ [Fe/H] <^ —3 in 
region A of Figure 2 indicates that the initial/prompt Fe 
production had diverse yields and/or was sufficiently ex- 
tended in time so that [Fe/H]~ —3 represents the sum of a 
number of individual events. The Ba abundances in region 
A could be attributed to production by the initial/prompt 
Fe source which would be small compared with that by 
a single H event. However, these data could also be ex- 
plained by a mixing scenario with no Ba production by the 
initial/prompt Fe source (see §4). The onset of the corre- 
lation between [Fe/H] and abundances of heavy r-elements 
can be estimated as follows. By attributing ^ 1/3 of the 
solar Fe inventory (Timmes, Woosley, & Weaver 1995) to 
the Type H supernovae associated with L events, we expect 
a single L event to result in logeL(Fe) 5.0 corresponding 

to [Fe/H] 2.5 if we take Al ~ 10^ yr (~ 100 L events 

are then responsible for the part of the solar Fe inventory 
contributed by Type II supernovae). Therefore, we ex- 
pect that a correlation between [Fe/H] and abundances of 
heavy r-elements would be established through addition 
of H and L events over a few 10^ yr during which the ISM 
was sufficiently enriched by L events to [Fe/H] ^ —2.5. 
Indeed, data in Figures 1 and 2 (region C) show that such 
a correlation exists at [Fe/H] ^ —2.5. 

4. DISCUSSION AND CONCLUSIONS 

There is a basic issue of what plausible mechanism 
could account for the initial or promptly-generated Fe in- 
ventory. The prompt Fe production in a gas clot must 
have lasted for only a narrow time interval (<C Ah ~ 
10^ yr) and then was greatly diminished. A possible 
mechanism would involve Fe production by supermas- 
sive stars with no significant coproduction of heavy r- 
nuclei. Let us consider a pristine gas clot of mass Mq 
from which two populations of stars could be made: su- 
permassive stars with very short lifetimes (^ 10^ yr) and 
less massive stars with lifetimes longer than ^ 10^ yr. 
So Mpit) + Mp{t) + Mg{t) = Mo, where Mpit), Mp{t), 
and Mg{t) represent the masses stored in supermassive 
("Fat") stars, less massive ("Petite") stars, and gas at 
time t. Taking the birth rates of both types of stars to 
be proportional to Mg{t) and assuming that supermas- 
sive stars were born at about the same rate as they were 
destroyed, we have Mpit) « KpMgit) - MF{t)/fF » 0, 
where fp is the average lifetime of supermassive stars, and 
Mp(t) = KpMg{t). Thus, Mp{t) « KpMofpexpi-Kpt) 
and Mp{t) « Mo[l - cxp{-Kpt)]. Therefore, if the birth 
rate for less massive stars was sufficiently high to deplete 
the gas over a few 10^ yr, then the population of super- 



massive stars would decline on the same timcscalc. This 
would provide some diversity in Fe abundances in stars 
formed at very early times but effectively truncate further 
addition of Fe from supermassive stars. 

Stars with masses 11 <J M/Mq <^ 40 are considered to 
become supernovae. Their lifetimes range from 6 x 10^ yr 
to 2 X 10^ yr (see Fig. 1 of Timmes et al. 1995). Meynet et 
al. (1994) gave a lifetime w 3 x 10^ yr for a 120 Mq star. 
Therefore, prompt Fe production in our model must be as- 
sociated with siipermassive stars of at least a few 100 Mq. 
These stars are assumed to produce no significant amount 
of heavy r-nuclei but sufficient Fe to give [Fe/H] ~ —3 
over a few 10^ yr. As speculated above, the termination 
of prompt Fe production was caused by rapid depletion of 
gas through storage in less massive stars over a few 10^ yr. 
Thus the initial star formation rate in a pristine gas clot 
of mass ~ 10^ Mq is required to be ~ IMq yr^^. When 
extrapolated to the whole Galaxy, this corresponds to an 
initial rate ~ 10''' Mq yr~^, much higher than the average 
value of ~ 10 Mq yr~^ over Galactic history. On a longer 
timescale, stars with 11 ^ M/Mq ^ 40 became super- 
novae and began to enrich the ISM. These supernovae are 
of two types: high frequency H events and low frequency 
L events. Within a standard mass of hydrogen, H and L 
events occur on timescales ~ lO'' yr and ~ 10^ yr, respec- 
tively. The H events produce heavy r-elements but no Fe 
and this resulted in a wide range of abundances for heavy 
r-elements (e.g., Ba and Eu) at [Fe/H] ^ —3. A correla- 
tion between [Fe/H] and abundances of heavy r-elements 
was established later through addition of H and L events 
over a few 10^ yr when Fe production by L events over- 
whelmed the prompt Fe inventory. The observed onset of 
this correlation at [Fe/H] ~ —2.5 is consistent with the 
expected Fe yield of L events. The L events also produce 
light r-elements such as Ag. Production of these elements 
by the hypothesized prompt Fe source is unknown. It was 
argued that the first stars of H-He composition would be 
very massive (e.g., Truran & Cameron 1971). These stars 
would provide elements heavier than He to the ISM (e.g., 
Ezer & Cameron 1971). How much Fe would be produced 
by these stars and how this Fe would be mixed with the 
ISM require more investigation. 

It is likely that supermassive stars would blow up their 
parent gas clots, thus preventing further star formation. 
In this case they would provide Fe to the general ISM. We 
speculate that once [Fe/H] ^ —3 was reached, supermas- 
sive stars could no longer be produced and less massive 
stars would form instead. In this scenario an initial Fe 
inventory could be provided without requiring a high star 
formation rate in the gas clots before the H and L events 
occurred. We note that the cut-off at [Fe/H] ~ —3 may 
correspond to a condition when sufficient amount of ele- 
ments heavier than He was provided by supermassive stars 
to permit adequate cooling of aggregating matter for less 
massive (^ 40 Mq) stars to form. A supermassive star of 
a few IOOA'/q must produce a few Mq of Fe in order to 
give [Fe/H] 3 to a clot of ~ 10^ Mq. 

Region A of Figure 2 deserves special attention. We 
are faced with Ba abundances below the production by a 
single H event at -4 <^ [Fe/H] <^ -3. These Ba abun- 
dances could be attributed to minor heavy r-element pro- 
duction by supermassive stars. However, supernova explo- 
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sions could drive gas outflows from a clot. The enriched 
gas could then mix with the pristine gas in other c;lots. In 
this way Ba abundances below the production of a single H 
event for [Fe/H] < —3 could be obtained. For example, the 
abundances in the star at the left corner of region A could 
be explained by a mixture of gas outflow after an H event 
with the gas in a pristine clot (mixing ratio ~ 1:10). We 
have avoided the complexities of mixing and exchange in 
our previous discussion based on a standard mixing mass. 
Subsequent models must address these issues. 

The chronometric interpretation of [Fe/H] is complex. 
Condensation of matter to form stars during the early 
evolution of the Galaxy is expected to have been greatly 
extended in space and time (~ 10^-10^ yr). This means 
that clots of baryonic matter formed within the Galaxy 
at widely disparate times. The model proposed here only 
requires the following sequence of events to occur within a 
clot: an initial Fe inventory from or prompt Fe enrichment 
by supermassive stars, enhancement in heavy r-elements 
by H events, and enrichment in Fe and light r-elcmcnts 
by L events. So long as different clots of baryonic matter 
within the Galaxy underwent the same evolution, this se- 
quence of events can be established independent of which 



stars observed today represent the same initial clot. How- 
ever, this same; sciqucncc of events might have started at 
widely different times within different clots. In this sense, 
until chemical enrichment became essentially uniform on 
the Galactic scale, stellar abundances at low metallicities 
would only reflect relative time in the above sequence. Fi- 
nally, if galaxies formed ^ 10^ yr after the Big Bang, the 
epoch of supermassive star formation discussed here would 
correspond to a redshift of at least 2; ~ 4 [for which the age 
of the universe is - 10^"(1 + z)"^/^ _ iq9 y^.]^ rj^^jg ^.p^^j^ 
is earlier than the one probed by recent abundance obser- 
vations at "high" redshifts {z ~ 3, Pettini et al. 1997). 
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Fig. 1.— Data (stars: Gratton & Sncdcn 1994; squares: McWilliam ct al. 1995; McWilliam 1998; triangles: Sneden et al. 1996; 1998) on 
[Ba/Eu]r = log(Ba/Eu) — log(Ba/Eu)0_r and loge(Eu) versus [Fe/H] for low-metallicity stars. The Ba/Eu ratio is near the solar r-process 
value over the range of [Fe/H] shown. But there is a wide dispersion in loge(Eu) at —3.1 < [Fe/H] < —2.5. This dispersion disappears for 
[Fe/H] > -2.5. 

Fig. 2. — Data (symbols as in Fig. 1) on logc(Ba) versus [Fc/H] for low-metallicity stars. The range of [Fe/H] extends about one dex below 
that for the existing data on Eu. There is a wide dispersion in loge(Ba) at —3.1 < [Fe/H] < —2.5. Three regions of abundance evolution 
are schematically shown: production by the initial/prompt Fe source (A), addition of high frequency non-Fe-producing H events (B), and 
mixture of H and low frequency Fe-producing L events (C). 
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